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A method is proposed for combining muhi-frequency maps in order to produce a catalogue of 
extragalactic point sources using data from future high-precision satellite experiments. We have 
found the optimal way for combining maps at different frequencies in order to maximize the signal 
(point sources) to noise (rest of the signal) ratio. Our approach is a natural multi-frequency gen- 
eralization of the band-pass filter introduced by Tegmark and de Oliveira-Costa (1998). We show 
that combination of different frequency maps gives us the possibility of creating a more complete 
catalogue of point sources. 
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I. INTRODUCTION 

In the next few years, the new generation of CMB 
experiments (MAP and PLANCK) will provide all-sky 
maps of the cosmic microwave background (CMB) with 
high resolutions and sensitivities. Theoretical investi- 
gations of different kinds of foregrounds, instrumental 
noise properties, scan strategies and related topics for 
the future PLANCK mission are now in progress. The 
PLANCK LFI and HFI instruments will be able to mea- 
sure the CMB anisotropy and polarization using 10 fre- 
quency channels that cover the frequency range 30 - 857 
GHz (Mandolesi et al.,1998; Maino et al.,1999). The 
maps provided by PLANCK will contain contributions 
from various physical components, including the primary 
and secondary CMB signals and foregrounds such as free- 
free and synchrotron emission, the Sunyaev-Zeldovich ef- 
fect, galactic dust and extragalactic point sources. One 
of the important goals of this experiment, apart from 
the separation of the CMB from the remaining parts of 
the signal, is the creation of a point source catalogue. 
The aim of this paper is to show how one may opti- 
mize the strategy for producing the point source cata- 
logue by using the multi-frequency properties of the ra- 
dio/submillimetre sky observations. 

The effect of point sources on satellite observations has 
been widely discussed in the literature (see the review by 
Vielva et al. (2001), and references therein). Many au- 
thors have developed various linear and non-linear meth- 
ods for point source extraction. Hobson et al. (1999) sug- 
gested the combination of the maximum- entropy method 
(MEM) with Mexican Hat wavelet filtering in order to 
separate all the physical components of the microwave 
ski including extragalactic point sources and using differ- 
ent frequency channels. Tegmark and de Oliveira-Costa 
(1998; hereafter TO-98) described the linear filtering of 
point sources in the form of a band-pass filter in order 
to render point sources on each frequency map more vis- 
ible relative to the background of other physical compo- 
nents. Cayon (2000), Sanz et al. (2000) and Vielva et al. 
(2000) applied the Mexican Hat wavelet method to de- 
noising the CMB maps and for production of the Planck 
point source catalogue. Naselsky et al. (2000), Chiang 
et al. (2001) proposed the approach of amplitude-phase 
analysis and applied it to simulated maps of the CMB 
with foregrounds, point sources and pixel noise included 
for the accurate determination of the locations of bright 
point sources and separation of these sources from the 
rest of the signal by an iteration technique. All of the 
methods mentioned above can complement each other in 
the framework of the future highly sensitive CMB ex- 
periments because of the different sensitivities of each 
method to different properties of the signal. 

In this paper we focus our attention on the combina- 
tion of maps with different frequencies in order to maxi- 
mize the signal {S- point sources) to noise (iV- rest of the 
signal) ratio. We show that co-adding of harmonics from 



different maps with specific weights leads to an increase 
in the accuracy of the point source subtraction. The aim 
of our paper is to produce the best possible combination 
of all frequency channels in order to construct a map in 
which the presence of point sources of some particular 
population is most visible. This procedure allows us to 
obtain a more complete and accurate catalogue of point 
sources than can be obtained by independently investi- 
gating each frequency channel (TO-98). 

The layout of our paper is as follows. In Section 2 we 
introduce the definition of the 'best' map combination, 
derive the analytical solution for the maximum possible 
signal-to-noise ratio and make some numerical predic- 
tions for the Planck experiment. We also present a pos- 
sible special combinations of the multifrequency maps. 
We briefiy summarize our results in Section 3. 

II. OPTIMAL COMBINATION OF DIFFERENT 
FREQUENCY CHANNELS 

In this section, we derive the optimal way for using the 

different frequency channels to produce the best possible 
map for detection and subsequent subtraction of extra- 
galactic point sources. The idea behind our approach is 
quite simple and transparent. We define the 'best' map as 
the combination of different frequency maps which max- 
imizes the signal-to-noise ratio. In our case, the signal is 
the signal from point sources. The noise consists of other 
physical components including pixel noise, CMB, dust, 
synchrotron and free-free emission. All of these compo- 
nents limit our capability to detect point sources. 

Optimization 

Suppose that we have N different maps of the same 
region of the sky made using N different channels with 
frequencies z/j, i = 1, .., N . The flux of each point source 
= S{v^) is a function of frequency and, therefore, 
differs from one channel to another, while the position 
of this source is, obviously, the same for all maps. In 
this subsection, we consider for simplicity one particular 
population of point sources with the same dependence of 
fixixes on frequency. The contribution of point sources 
to the the temperature on the CMB map in the j-th fre- 
rjuency channel can be written as a set of Dirac delta 
functions with coefficients proportional to the fiux of 
the source at position fj: 

^Ti,=Y,si5{n,T) (1) 

i=l 

Here, Np^ is the total number of point sources. The tem- 
perature fluctuations on the map convolved with the an- 
tenna beam are therefore: 

AT^" (f ) = E slB^ (fl - f ) + E aimYi^ (f ) = 

Im 
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FIG. 1. Signal-to-noise ratio (in arbitrary units) for each 
channel shown separately after applying the TO-98 filter 
(dashed lines) and for the combination of 10 Planck maps 
(solid line). The calculation is performed for the population 
of radio sources with spectral parameter a taken from Table 
1. The dotted line shows the same ratio, but for a break in 
the spectrum between 143 and 217 GHz. 
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FIG. 2. Same as fig.l but for the population of FIR sources 

^ 7 = Afps(fi)/CTtot, 



where (f ) is the beam for the j-th channel, the am- 
plitude of the point source si is proportional to the flux 
and depends on the antenna width: ^ S^/Q"^. The 
coefficients aj^ have zero mean and variance {{0,1^)"^) = 
^imiHm)^ + ^pix (pi^^l noise is not convolved with the 
beam). The set of harmonics T^-^, j = 1,..,N can be 
considered as the N-dimcnsional vector: Tim- The aim 
of our paper is to combine N given maps with known 
values of the expected noise to construct the map with 
the best signal-to-noise ratio. In order to produce such a 
map, we introduce the vector-filter /j^ with components 
■firm .., and consider the combination of N maps in 
the following form: 

AT = E tmTlmYlm (f) = E TlmYlm (f) (3) 
Im Im 

~ N . . 

where T/^ = E fim'^lm- consider a symmetric an- 

i=i 

tenna beam and subsequently drop the subscript m be- 
cause {{Tim) ) = {{Ti) ). Thus the filter that we use 
should depend on 1 only: /; = /;^. Our approach can 
readily be generalized to the case of an anisotropic beam, 
but this is beyond the scope of the present paper. Our 
problem now is to find vectors /; for each I so that the 
ratio 7 = signal /noise is maximized: 



Here /?; is the vector with components sjV , j = 

1, ..A^ (6; is the beam) and Mi = [mf^'-'^] is the covariance 

Ji h \ 



(5) 



where (5j^ is the Kronecker symbol. The ratio 7 from 
equation (4) is maximized if we choose the filter to be: 



(6) 



Consequently, one can conclude that point sources with 
some particular dependency of flux on frequency are most 
clearly seen in the map that is the combination of maps 
taken with these specific weights of harmonics. If only 
one channel is under consideration, then this formula be- 
comes especially simple and is exactly the same as the one 
obtained by Tegmark and de Oliveira-Costa (TO-98): 



fi 



bi Ci + Cpix 



(7) 



Using equations (4) and (6) , one can get the final result 
for the maximum possible signal-to-noise ratio 7, which 
is achieved in the combined map: 
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1) and increasing to — oo at higher frequency. For the 
far-IR population, this spectral index has the opposite 
behavior and we choose a to be « 2.7 for 857-352 GHz 
with practically no contribution from these sources to the 
signal at lower frequencies. Of course, our estimates are 
applicable only for these two specific populations. Each 
particular population of point sources should be inves- 
tigated separately. This means that we should perform 
the combination of channels with specific filters /; which 
correspond to the appropriate spectral index a for this 
population. 

Using the results of the previous subsection, one can 
calculate the signal-to-noise ratio for each channel sepa- 
rately 7(j) (if only the j-th channel is under consideration 
(TO-98)) and for the combination of N channels. From 
eq. [7,8] one finds: 



7' = Em + ihf. 



(10) 



FIG. 3. Signal-to-noise ratio for 143 GHz and 217 GHz 
channels after performing TO-98 filtering and for the com- 
bination of channels. Lower straight lines show signal/noise 
for initial maps (in units of atot) while upper straight lines 
show this ratio after TO-98 filtering. Dotted fines: 217 GHz, 
dashed lines: 143 GHz, solid fine: combination of maps. 



^(2Z + l)/3fMr^A 



1/2 



(8) 



The contribution to signal /noise in each / is obviously 
higher for the combination of maps than for each map 
separately and consequently 7 > 7(j). In fig. (1,2) 
we show this ratio for two populations (radio and FIR) 
as described above. A significant improvement can be 
achieved by considering both populations, 7 = A x 
max[y{j)], j=l,..,10. For radio source, s we found A « 4.4 
and for FIR A is about 1.2. Therefore, a more complete 
catalogue can be produced for this population of sources. 

Special cases and a possible combinations 



Estimates for the Planck experiment 

A multi-frequency analysis gives us the possibility of 
detecting point sources in coincidence with different fre- 
quency maps. In our calculations, we assume the fre- 
quency dependence of the intensity for some particu- 
lar point source population in the form of a simple 
power law: p". Multiplying this intensity by the conver- 
sion factor between surface brightness and temperature 
(Tegmark and Efstathiou, 1996) and taking into account 
the resolution 8, one can obtain the dependence of the 
point source temperature on frequency: 

^^,^^^^ s^nH.,/,6.SGH^) ^ ,- = l,..,10. (9) 

where j denotes the Planck channel number. Channels 
are numbered in order of increasing frequency (see Table 
1,2). Two different populations have been chosen: radio 
and far-IR sources. Up to about 200 GHz, we can ex- 
pect mostly radio sources to be detected, while at higher 
frcquenies far-IR sources dominate. For the radio source 
population, we took the spectral index a from (Vielva et 
al., 2001) to be close to zero up to 100-200 GHz (table 



Let us for simplicity consider a special case by dealing 
with only two channels with different frequencies: !^i,f2- 
Since the spectral behavior of different populations of 
point sources is not well known, we introduce the pa- 
rameter p. This parameter contains information about 
amplitudes of temperature fluctuations due to the same 
point sources in different channels: 

ATpi, = pal,, ATI, = (1 - P)^tot (11) 

In fig. 3 we show the signal-to-noise ratio that can be 
found applying TO-98 filters for the 143 and 217 GHz 
channels separately and for the combination of these two 
maps with different possible ratios of amplitudes AT^, 
and ATpg. We choose these two channels because sig- 
nificant steepening or even a break in the point source 
spectrum due to synchrotron emission is expected over 
this range of frequencies, so it is hard to predict the ex- 
act value of p. 

For smaller values of u quite a large number of point 
sources are expected to be found in coincidence for dif- 
ferent channels. It is interesting to investigate how does 
our filter look like if we combine for example the 70 GHz 
and 100 GHz (LFI) maps. 
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This means that our filter removes contributions from 
the CMB and other physical components at scales larger 
that the typical scale of the point sources. 

For I values that correspond to scales comparable with 
the antenna beam, the pixel noise is dominant and: 
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(15) 



Therefore the components of the filter /; arc close to the 
filters found in TO-98 for each channel separately (see 
fig.4): 



pix 



pix 



(16) 



These components are both positive and add to form 
peaks that correspond to point sources in the combined 
map. 

III. CONCLUSIONS. 



FIG. 4. Components of the filter /( for combining the 
best possible map from the 70 and 100 (LFI) GHz channels 

(solid lines) in the comparison with TO-98 filters for each 
map taken separately (dashed lines). The filter is given in 
arbitrary units. 

According to formula (6) one can obtain the filter /; 
with components f^, ff: 



Ji — 



det{Mi) 
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det{Mi) 



(12) 



where det{Mi) = ml^mf^ — (m^)^. Using eq [5] we can 
investigate the behaviour of these two components at dif- 
ferent scales. 

For small Ps that correspond to scales much larger 
than the antenna resolution, the instrumental noise is 
negligible (at least for the low frequency channels: v < 
353GHz) and: 



C? 



(13) 



m 



12 



Substituting [13] in [12] we obtain: 



Wc have presented an analysis of microwave sky obser- 
vations based on the combination of channels with dif- 
ferent frequencies in order to separate extragalactic point 
sources for subsequent production of a point source cata- 
logue. We made some numerical estimates for the future 
Planck experiment and showed that a multi-frequency 
analysis allows us to increase the signal-to-noise ratio by 
a factor of 1.2-4.4 to create a more complete and precise 
catalog of point sources. 



fl - deL{M) - f^lV^) ' 

AW) • 



det{M) 



(14) 



Therefore, harmonics with small / from different channels 
are taken with the opposite sign. Thus the contributions 
from the first and second maps to the combined map 
compensate each other. Finally the total signal in the 
combined map turns to be close to zero at large scales. 
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[^KjtlZ ) 


[a) 


(initial) 


i *j-yo ) 


(combination) 


1 


30 


0.08 


1.000 


3.083 


3.066 


2 


44 


0.08 


0.986 


3.091 


5.894 


3 


70 


-0.10 


0.935 


2.495 


8.579 


4 


100(L1< 1) 


-0.16 


0.947 


2.490 


10.862 


5 


lOO(HFI) 




0.872 


3.622 


13.026 


6 


143 


-0.55 


0.766 


3.429 


14.858 


7 


217 


-0.55 


0.952 


3.157 


15.747 


8 


353 


-0.55 


0.895 


1.886 


15.834 


9 


545 


— oo 


0.000 


0.000 


15.847 


10 


857 




0.000 


0.000 


15.925 


TABLE I. 


Results of combination of 10 Planck channels for the population of radio sources with given spectral behavior. The 


frequency of each channel is indicated in column 2 for both LFI and HFI. The 3-rd column 


gives spectral indices as calculated 


by Viclva ot al. (2001). — oc denotes 


spectral break. Fourth column shows the initial signal-to-noisc ratio in units of a^^^ under 


the condition that the point source 


in the 1-st channel has la"tot 


level. The fifth column 


gives this ratio after applying the 


TO-98 filter to each channel separately. In the last column we show this ratio for the best combination of channels from 1 to j. 


This means 


that the 2-nd row gives the result for the combination of the 1-st and 2-nd channels, 3-d 


means the combination 


of 1,2,3 channels and so on. 










Channel 


Frequency 


spectral index 


signal/noise si 


gnal/noise 


signal/noise 


number (j) 


(GHz) 


(a) 


(initial) 


(TO-98) 


(combination) 


10 


857 


2.46 


1.000 


6.281 


6.281 


9 


545 


2.71 


0.811 


2.257 


6.575 


8 


353 


2.71 


0.551 


1.162 


6.960 


7 


217 


oo 


0.000 


0.000 


7.405 


6 


143 




0.000 


0.000 


7.460 


5 


lOO(HFI) 




0.000 


0.000 


7.468 


4 


lOO(LFI) 




0.000 


0.000 


7.470 


3 


70 




0.000 


0.000 


7.470 


2 


44 




0.000 


0.000 


7.470 


1 


30 




0.000 


0.000 


7.470 



TABLE II. Same as Table 1 but for FIR sources population. The channels arc shown in the opposite order. Initial signal/noise 
ratios are given under the assumption that the amplitude of the source in the 10-th channel is lal^t- 
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